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Abstract. The global distribution of dust column burden de- 
rived from MODIS Deep Blue aerosol products is compared 
to NH 3 column burden retrieved from IASI infrared spectra. 
We found similarities in their spatial distributions, in partic- 
ular their hot spots are often collocated over croplands and 
to a lesser extent pastures. Globally, we found 22 % of dust 
burden collocated with NH 3 , with only 1 % difference be- 
tween land-use databases. This confirms the importance of 
anthropogenic dust from agriculture. Regionally, the Indian 
subcontinent has the highest amount of dust mixed with NH 3 
(26%), mostly over cropland and during the pre-monsoon 
season. North Africa represents 50 % of total dust burden but 
accounts for only 4 % of mixed dust, which is found over 
croplands and pastures in Sahel and the coastal region of the 
Mediterranean. In order to evaluate the radiative effect of this 
mixing on dust optical properties, we derive the mass extinc- 
tion efficiency for various mixtures of dust and NH 3 , using 
AERONET sunphotometers data. We found that for dusty 
days the coarse mode mass extinction efficiency decreases 
from 0.62 to 0.48 m 2 g -1 as NH 3 burden increases from 0 
to 40mgm~ 2 . The fine mode extinction efficiency, ranging 
from 4 to 16 m 2 g _1 , does not appear to depend on NH 3 con- 
centration or relative humidity but rather on mineralogical 
composition and mixing with other aerosols. Our results im- 
ply that a significant amount of dust is already mixed with 
ammonium salt before its long range transport. This in turn 
will affect dust lifetime, and its interactions with radiation 
and cloud properties. 


1 Introduction 

Mineral dust affects climate by absorbing and scattering so- 
lar and terrestrial radiation, by modifying cloud properties, 
and by interacting with ozone chemistry and the biosphere 
(Forster et al., 2007). The amplitude of these effects depends 
on the size, optical and chemical properties of dust. These 
properties may be considerably modified when dust is mixed 
with other aerosol types. In the case of dust mixed with am- 
monium salts, ( 1 ) the amount of solar absorption will de- 
crease (Bauer et al., 2007), (2) the activation of cloud con- 
densation nuclei will increase (Levin et al., 1996), (3) the ac- 
tivation of ice nucleation sites will decrease (Sullivan et al., 
2010 ), (4) the recycling of NO x is suppressed reducing pol- 
lution influence on surface ozone (Fairlie et al., 2010), and 
(5) iron dissolution which enhances ocean fertilization will 
increase (Meskhidze et al., 2003). Therefore, it is crucial to 
determine the mixing state of dust in order to better under- 
stand its effects on climate. 

In situ observations have shown that ammonium salts com- 
monly accumulate on dust particles (Galy-Lacaux et al., 
2001; Jordan et al., 2003; Sullivan et al., 2007; Shi et al., 
2008). The coating of dust particles by sulfate occurred by 
in-cloud processes, by which reactions with sulfuric acid 
(H 2 SO 4 ) initiate the deposition of sulfate on the surface 
of the particles (Levin et al., 1996). Other mixing pro- 
cesses include coagulation of sulfate particles or oxidation 
of sulfur dioxide (SO 2 ) to sulfate on the particle surface 
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in the presence of ammonia (NH3). The latter process in- 
volves titration of NH3 by H2SO4 to form ammonium sulfate 
(NH4)2S04, ammonium bisulfate (NTLOSC^, or triammo- 
nium hydrogen disulfate (NH4)3H(S04)2. Similarly, in pres- 
ence of nitric acid (HNO3), NH3 will produce ammonium 
nitrate (NH3NO3) (Usher et al., 2003). Sulfate coatings on 
dust undergo efflorescence (crystallization) at different rela- 
tive humidities depending on the size and mineralogy of dust 
(Usher et al., 2003). 

By analyzing Asian dust, Sullivan et al. (2007) showed 
that transport pathways of dust through source of pollu- 
tants is more important than reaction kinetics in determining 
which species accumulate on dust. The scarcity of in-situ ob- 
servations makes it difficult to extrapolate the predominance 
of dust mixing with pollutants at global scale, such that the 
impact of such mixing on climate is evaluated using poorly 
constrained models (Lesins et al., 2002; Bauer et al., 2007). 

Using the Moderate Resolution Imaging Spectroradiome- 
ter (MODIS) Deep Blue satellite products, Ginoux et al. 
(2010) have retrieved global dust sources at 0.1° resolution. 
By comparing their location with land-use fraction, Ginoux 
et al. (2012) have calculated that 25 % of dust is emitted from 
agriculture, and called such dust as anthropogenic. Mean- 
while, Clarisse et al. (2009) have retrieved the global dis- 
tribution of ammonia (NH3) from the Infrared Atmospheric 
Sounding Interferometer (IASI) infrared spectra and have 
shown that hot spots are often associated with agriculture. 
There is a remarkable similarity in the observed distribution 
of these hot spots and anthropogenic dust sources, which may 
provide the first evidence of co-located pollutants and dust 
over their sources (see Sect. 3). This would first confirm the 
existence of anthropogenic dust (Sect. 4), once our results 
have been validated (Sect. 5), and second that anthropogenic 
dust has different chemical and optical properties (Sect. 6). 
The implications of these results on our understanding of the 
effects of dust on climate are discussed in Sect. 7. 


2 Satellite data 

2.1 Dust column burden 

The MODIS instrument on the polar orbiting Aqua satellite 
measures daily backscattering radiances globally since 2002. 
The MODIS Deep Blue algorithm employs radiances from 
the blue channels of MODIS instruments. At these wave- 
lengths the surface reflectance is very low so that the pres- 
ence of aerosol is detected by an increase of total reflectance 
and enhanced spectral contrast (Hsu et al., 2004, 2006). The 
values of AOD and single scattering albedo (< co ) at 412, 470, 
550 and 670 nm, and the Angstrom exponent (a) between 
412 and 470 nm are retrieved at the pixel level over bright 
surfaces (reflectance at 550 nm greater than 0.15), and then 
averaged on a 10 by 10 km grid. The data are aggregated into 
granules forming the Level 2 data. These data are interpo- 


lated by surface-weighted averaging on a 0. 1 0 by 0. 1 0 regular 
latitude/longitude grid. Ginoux et al. (2012) extracted Dust 
Optical Depth (DOD) from MODIS collection 5.1 Level 2 
(M-DB2) aerosol products by requiring the values of the 
Angstrom exponent to be less than 0, and the single scat- 
tering albedo to increase with wavelength. By analyzing 7 yr 
of daily global DOD, Ginoux et al. (2012) also found that 
globally 24 % of dust is emitted from sources associated with 
land-use fraction greater than 30 %, which they defined as an- 
thropogenic. They also found very large variations between 
continents and seasons. 

They validated M-DB2 DOD by comparing with 
AERONET data. They showed significant correlation be- 
tween M-DB2 and AERONET data, with root mean square 
difference of 0.26 and mean absolute difference of 0.24. 
The largest biases were observed in California and Australia 
where M-DB2 largely overestimates DOD. 

In the present study, we use daily global M-DB2 DOD 
from April 2009 to March 2010 and convert it into 
column burden. The conversion is performed following 
Ginoux et al. (2001) who expressed DOD (r) and mass bur- 
den M by the following relation, 



P r eff 

T 

(2 ext 


1 

-r 


( 1 ) 


where r e ff = 1.2 x 10 -6 m is the effective radius, p = 
2600 kgm -3 is the density of dust, Q ext = 2.5 is the extinc- 
tion efficiency at 550 nm and foral.2xl0 -6 particle radius, 
e = 0.6 m 2 g -1 is the mass extinction efficiency, and r is dust 
optical depth at 550 nm. The selection of r e ff is based on the 
effective radius of the total size distribution measured at Sal 
Island (Cape Verde) and Banizoumbou (Niger) by Tanre et al. 
(2001). As the aerosol optical depth increased, they observed 
that r e ff reaches 1 .2 x 10 -6 m at these two sites. 

Such relation was also used by Kaufman et al. (2005) 
to convert DOD to M. Using r e ff = 1.7 x 10 -6 m, 
p = 2600 kgm -3 and Q ex t = 2.2, they obtained 
€ = 0.37 =b 0.05 m 2 g _1 , which is 40% lower than our 
value. On the other hand, our value is closer to the me- 
dian value (0.72 m 2 g -1 ) derived from global dust models 
(Huneeus et al., 2011). The range of values among these 
models is quite large, varying from 0.25 to 1.28m 2 g _1 . 
Similar range has been obtained from theoretical calcula- 
tions (Hansell et al., 2011), with higher values associated 
with non-spherical particles. Also, values derived from field 
campaigns vary from 0.6 d= 0.1 at Barbados (Li et al., 1996) 
to 1.09 =b 0.4 at Cape Verde (Chen et al., 2011). There is 
clearly a need to better constrain e, particularly over source 
areas where there are fewer data. This will be discussed in 
Sect. 5. 
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Fig. 1. Global (lower panel) distribution of M-DB2 dust column burden (gm 2 ; grey shading) and IASI NH3 column burden (mgm 2 ; color 
shading), on 16 May 2009. The upper panels provide a zoom over 6 regions. 


2.2 NH 3 column burden 

Infrared radiances have been measured daily since the end of 
2006 by the IASI instrument on board the EUMETSAT’s po- 
lar orbiter MetOp. Clarisse et al. (2009) calculated the bright- 
ness temperature difference between radiances and converted 
it to column burden of NH 3 . The first global distributions 
were obtained by averaging one year of measurements and 
28 hotspot regions were identified. Here we use an improved 


product only available from April 2009 to March 2010. For 
this product NH 3 column burden is retrieved using a full Op- 
timal Estimation scheme, as in Clarisse et al. (2010a), and 
similar to the processing algorithm used for other trace gases 
(Hurtmans et al., 2012). The retrieval is only carried out for 
the observations with a detectable spectral signal of NH 3 . 
Hence, very low concentrations will not be measured. Simi- 
larly, when the contrast in temperatures between the ground 
and the atmosphere is weak, the IASI sensitivity to surface 
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concentrations is close to zero, and even significant NH 3 
concentrations might be missed. While the IASI optimal es- 
timation retrievals account for emissivity features by using 
temporal and spectrally resolved emissivity data (Zhou et ah, 
2011 ), some residual contamination due to emissivity fea- 
tures of windblown dust (Clarisse et ah, 2010a) cannot be 
excluded. To minimize these, only those NH 3 retrievals with 
root mean square of the observed minus fitted spectrum be- 
low 2.8 Won -1 sr -1 were retained. This threshold value is 
close to the IASI instrumental noise in the relevant spectral 
range (Clerbaux et al., 2009). 

Here we use the daily NH 3 column burden retrieved from 
April 2009 to March 2010 on each pixel (12 km diameter at 
nadir). Monthly averages are computed on a 0.25° by 0.25° 
grid. The value in each grid cell was obtained by averaging 
all observations within that cell. In grid cells where no mea- 
surements were made a value of zero was assigned. Although 
some biases cannot be avoided, this approach minimizes ar- 
tifacts that would be introduced from interpolation of sparse 
data. 


3 Global distribution of anthropogenic dust and NH 3 

Figure 1 shows the column burden of M-DB2 dust, and IASI 
NH 3 on 16 May 2009. On that day significant amount of dust 
and NH 3 was observed to co-exist in different areas in the 
northern hemisphere. The main hot spots are in the south- 
ern part of the San Joaquin valley in California, the coastal 
region of Senegal, the Niger River Basin around Niamey 
(Niger), the Jordan River Basin, the Ganges River Basin, 
and the Horquin Sandy land in East China. All of these hot- 
spots have been identified as anthropogenic dust sources by 
Ginoux et al. (2012). This snaphsot of dust and NH 3 ex- 
amplifies their widespread co-existence over anthropogenic 
sources. The colocation of their sources suggests their mix- 
ing. However, the sources within the Niger River Basin are 
downwind of the Bodele depression, one of the most active 
dust source on Earth (Prospero et al., 2002). MODIS swath 
on Aqua did not cover the depression on that day, which ex- 
plains why Figure 1 does not show the Niger River Basin 
plume emanating from the Bodele on that day. Similar occur- 
rences of dust transport in elevated layers over Sahel, while 
NH 3 resides in the boundary layer, may falsely be interpreted 
as dust mixedintroduce with NH 3 . On the other hand, it has 
been suggested from modeling studies with support from air- 
craft data that NH 3 may, in some cases, be transported in the 
free troposphere (Wang et al., 2008). Using vertical profile 
of aerosols from lidar in space satellite instruments may par- 
tially help, but their narrow field of view makes seldom their 
collocation with location specific small dust events. 

We would like now to estimate the persistence of collo- 
cated sources an annual basis, and estimate the fraction of 
dust burden mixed with NH 3 for different types of land use. 


Klein Goldewijk (2001) (KG) and Ramankutty and Foley 
(1999) (RF) have derived the global distribution of the frac- 
tion of pasture and cropland per 0.1 and 0.5 degree gridcell, 
respectively. KH derived land-use type from statistical in- 
formation on national and sub-national level and population 
density, while RF used remote sensing ans statistical anal- 
ysis. As in Ginoux et al. (2010), we use here KG land-use 
database, but will also compare our results using RF data. 

We distinguish between natural and anthropogenic (agri- 
cultural) dust by assuming that the latter corresponds to dust 
over areas with more than 30 % land use. Natural dust is as- 
sociated with less than 30 % land use. We further separate 
between cropland and pasture dust depending on which one 
has the largest fraction within a grid cell. 

Figure 2 shows the annual mean column burden of M-DB2 
dust, and IASI NH 3 . We see that most of the “dust belt”, 
defined by Prospero et al. (2002), is covered with at least 
0.1 gm “ 2 of dust. The Sahara desert is essentially character- 
ized by natural dust while most of Sahel is associated with 
pasture except in Senegal, Burkina Faso, Niger and Nige- 
ria where cropland dominates. In the Middle East, cropland 
is confined to some areas of Mesopotamia. In Central Asia, 
most anthropogenic dust is associated with pasture, while 
dust over cropland covers most of India and East China. In 
the Southern Hemisphere, dust over cropland is found in Ar- 
gentina and Southeast Australia, and dust over pasture cov- 
ers Patagonia and Central and Northeast Australia. Over the 
United States, dust over cropland is confined to the southern 
part of the San Joaquin Valley of California and the Midwest, 
and dust over pasture dominates the Southwest. Major differ- 
ences of land-use fraction between KG and RF exist over the 
Saudi Arabian Peninsula and Australia, where land use is less 
than 30 % in case of RF database. 

Inventory of global ammonia emissions indicates that 
60% comes from agriculture, with 1/3 from fertilizer and 
2/3 from livestock sources (Beusen et al., 2008). Lifetime of 
gaseous NH 3 is about 22 h (Adams et al., 1999), which pro- 
duces well defined hot spots centered over its sources. This 
is clearly apparent in Fig. 2 where NH 3 distribution is cover- 
ing most cropland, which happens to be also source of dust 
(i.e. US Great Plains, Argentina, Sahel, Mediterranean coast, 
Mesopotamia, Uzbekistan, India, East China). The most pro- 
nounced feature is over the Indo-Gangetic basin, where dust 
and NH 3 column burden are greater than 0.5 gm~ 2 , and 
0.5 mgm -2 , respectively. 

The global annual mean dust burden from M-DB2 is 
4.2 Tg, which is four times lower than the median value of 
16 Tg derived by Huneeus et al. (201 1) from 14 models, with 
values ranging from 6.8 to 29.5 Tg depending on the model. 
There are two major reasons for this discrepancy: one is due 
to uncertainty in € in Eq. (1) and the other comes from M- 
DB2 retrieval over bright surface only. By excluding dark 
surfaces (i.e. oceans) transported dust is not included in our 
results. As we are interested to study dust over its sources, 
which are located over land, this is not an issue. Additionally, 
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Fig. 2. Global (lower panel) annual mean distribution of M-DB2 dust column burden (gm“ 2 ) over areas with less than 30 % (Natural Dust; 
red shading) and more than 30 % land-use -with an additional distinction between the prevalence of pasture over cropland (pasture dust; green 
shading) and cropland over pasture (cropland dust; blue shading), 0.5 (yellow contour) and 2.5 (orange contour) column burden (mgm -2 ) 
of IASI NH 3 , with the numbers 1 to 10 indicating the position of the AERONET sites of Fig. 4, and the upper panels providing a zoom over 
6 regions. 


we consider only one year while interannual variability of 
dust burden is of the order of 20 % (Ginoux et al., 2004). 


4 Anthropogenic dust mixed with NH 3 

In this section, we quantify the mixing of dust with NH 3 over 
agriculture at the global and continental scales. 


The percentage of anthropogenic dust mixed with NH 3 
over each landcover types is calculated by summing areas 
of each type where M-DB2 dust is mixed with detectable 
amount of IASI NH 3 (mass burden greater than 1 mgm -2 ), 
and dividing by the total surface of that cover type where 
M-DB2 dust is present. We consider both KG and RF land- 
use databases. We see in Fig. 3 that, for both databases, the 
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Table 1. Annual mean total burden of M-DB2 dust (M^u) i n units of ktons as well as the fraction (%) of dust mixed with NH 3 over all land 
cover types (/ tot ), cropland (/ crop ), pasture (/ past ), and with less than 30 % land-use (/ nat ), for the Klein Goldewijk (2001) and Ramankutty 
and Foley (1999) (in italics) landuse datasets, over 9 regions which boundaries are defined by their minimum and maximum latitude (Lat) 
and longitude (Lon). Percentage fraction greater than 20 % are indicated in bold font. 


Region 

Lon Range 

Lat Range 

M tot 

m du 



M nat 

m du 

ftOt 

ycrop 

ypast 

/ nat 

Global 

130° W-I6O 0 E 

60° S-60°N 

4184 

532 -579 

1622-7422 

2029-2183 

1-7 

22-21 

1-8 

4-4 

North America 

125° W-70° W 

20° N-50° N 

128 

9-15 

85-57 

33 -26 

11-77 

19-75 

11-77 

10-76 

South America 

85° W-60° W 

55° S-0° N 

70 

3-5 

28 -26 

38 -39 

6-6 

26-75 

3-77 

6-76 

North Africa 

20° W-35° E 

5° N-40° N 

1934 

141-755 

661-627 

1126-7742 

4-4 

11-76 

7-6 

2-2 

South Africa 

5° E-50° E 

35° S-5° N 

95 

2-7 

10-65 

24 -29 

5-5 

13-6 

5-5 

4-5 

West Asia 

35° E-60° E 

5° S-50° N 

1000 

12-70 

431 -343 

498 -587 

4-4 

8-7 

5-6 

4-3 

Central Asia 

60° E-100° E 

5° N-30° N 

334 

223 -248 

0-7 

111 -85 

26 -26 

34-32 

0-23 

12-76 

East Asia 

60° E-140° E 

30° N-50° N 

481 

63 -59 

222-206 

196-276 

14-74 

25-26 

15-76 

9-5 

Australia 

110° E-155°E 

45° S-10° S 

142 

6-3 

133 -78 

3-67 

1-7 

3-3 

1-7 

0-6 


Fraction of dust mixed with NH 3 over various landuse types 



Fig. 3. Annual mean fraction (%) of M-DB2 dust column burden 
mixed with IASI NH 3 over all land surfaces (blue), areas with less 
than 30% land-use (brown), cropland (red), and pasture (green) 
globally and over 8 continental regions, considering the Klein Gold- 
ewijk (2001) (KG) and Ramankutty and Foley (1999) (RF) land-use 
datasets, with the latter shown using darker colors. 


largest percentage of mixed dust with NH3 is over cropland 
globally, and for each sub-continents. 

Table 1 provides more details on the annual burden and 
fraction of dust mixed with NH3 over various continental re- 
gions and for KG and RF land-use databases. Globally, 7 % 
of dust is mixed with NH3, which represents 3 10 ktons of 
dust. If we consider only croplands, we found 21-22% of 
dust is mixed with NH3, which corresponds to 1 15-122 ktons 
or 3 % of global dust burden (4.2 Tg), depending on the 
land-use database. The largest fraction of dust mixed with 
NH3 (26 %) is observed over croplands of South and Cen- 
tral Asia subcontinent, with 32-34 % corresponding to 66- 
76 ktons, depending on the land-use database. The other area 
with significant amount of mixed dust is East Asia (67 ktons 
or 14% of global burden) with 25-26% from cropland. In 
North America, mixed dust represents only 14 ktons (11 % 
of North American dust burden) mostly from cropland (18- 
19%). North Africa represents 46% of global dust burden, 


but only 77 ktons (less than 2% of global burden) is mixed 
with NH3 with most of it over cropland. 

5 Comparison with AERONET data 

In order to evaluate M-DB2 derived column burden of 
dust, its value is compared with retrieval from inversion of 
Aerosol Robotic Network (AERONET) sunphotometer data. 
AERONET is a federated worldwide network of sunpho- 
tometers that are monitored and maintained at the NASA 
Goddard Space Flight Center (Holben et al., 1998). Dust par- 
ticles are essentially greater than 1 pm radius (coarse mode), 
but about 10 % of dust mass, mostly clay, is submicron (fine 
mode) (Ginoux et al., 2001). The parameters of the size dis- 
tribution for fine and coarse modes are retrieved by inversion 
of Almucantar measurements (Dubovik and King, 2000). 
Here we use the daily quality assured (Level 2) volume con- 
centration for the coarse mode from 1 April 2009 to the end 
of March 2010 in order to coincide with IASI data. 

The Angstrom wavelength exponent a is often used as an 
indicator of the dominance of one aerosol mode over the 
other. For example, using AERONET data, Eck et al. (1999) 
have shown that a values are generally lower than 0.5 in 
dusty environments and higher in polluted regions. To make 
sure that dust is the major aerosol contributing to the mea- 
sured extinction, we extract coarse mode volume concentra- 
tion when a < 0.5. The volume concentration is then mul- 
tiplied by the uncoated dry dust density (2600 kg m -3 ) to 
convert it into column burden. Samples collected in differ- 
ent areas of North Africa contain between 2 and 6 % sulfate 
in volume (Wagner et al., 2012), which suggest an error of 
the order of 10 % in the density. 

The AERONET retrieval algorithm in case of dust is de- 
scribed in details by Dubovik et al. (2006). Here, briefly, 
the algorithm assumes spherical particles, except when the 
Angstrom exponent is lower than 0.5, in which case spheroid 
models are used. For Angstrom less than 0.5, the estimated 
error is between 10 and 20 % for Level 2 data. 
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Fig. 4. Daily (small dots) and monthly (large dots) AERONET coarse mode column burden (black), M-DB2 dust burden (brown) in units 
of gm -2 (left vertical axes), and IASI NH 3 column burden (green) in units ofmgm -2 (right vertical axes) at 10 AERONET sites. The root 
mean square difference (RMSD) in units of gm -2 between monthly M-DB2 and AERONET mass column burden is provided in the upper 
right. 


Ten AERONET sites have been selected based on the di- 
versity of their geographical location over agriculture dom- 
inated areas. Satellite data within 25 km of the AERONET 
sites are used for this comparison. 

Figure 4 shows the time series of daily and monthly dust 
column burden from M-DB2 and AERONET and monthly 
IASI NH 3 column burden at the selected sites from April 
2009 to March 2010. The largest root mean square differ- 
ence (RMSD) between AERONET and M-DB2 dust burden 


appears in Niger at Banizoumbou (RMSD = 0.44 gm -2 ) and 
Zinder (RMSD = 0.41 gm -2 ), but the mean values are also 
the largest at these sites. The highest monthly dust burden 
(2.3 gm -2 ) among all ten sites is observed in Banizoum- 
bou (Niger) in March. African sites are characterized by two 
maxima of dust burden, one in March and the other in boreal 
summer (June- July- August). The seasonal variation of dust 
burden at the sites located in Sahel (south of 20° N) can be 
explained by the movement of the InterTropical Convergence 
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Zone (ITCZ) occupying its southernmost position (5° N) in 
December and its northernmost position (20° N) in July. As 
the ITCZ moves over these sites, dry air mass from the north 
is replaced by monsoon flow from the south bringing precipi- 
tation, which suppresses dust emission and removes dust par- 
ticles from the atmosphere. Sharp changes of wind direction 
and relative humidity have been observed in mid-October 
and early May at Banizoumbou, as well as a maximum of 
dust optical depth in March and June (Rajot et al., 2008). 
Summer maximum of dust load at the Kuwait site has been 
associated with increased surface winds in summer (Sabbah, 
2010). The RMSD (0.22 gm“ 2 ) is relatively low at that site, 
but is reaching 0.3 gm ~ 2 in Sede Boker (Israel). The summer 
maximum in Sede Boker has been associated with dust orig- 
inating from the Arabian Peninsula and the Sahara (Israele- 
vich et al., 2003). At the two sites in the Indo-Gangetic basin 
(Karachi and Jaipur), there is a pronounced maximum of dust 
load during the pre-monsoon season with monthly dust bur- 
den reaching up to 1.5 gm~ 2 . This seasonal variation cor- 
responds to the climatological (2001-2009) cycle of coarse 
mode optical depth measured with sunphotometer over the 
Indo-Gangetic Plain (Giles et al., 201 1). In China, dust activ- 
ity peaks in spring (Prospero et al., 2002) but M-DB2 value 
underestimates AERONET data by a factor 2 in March. This 
large discrepancy may be due to the low number of observa- 
tions during that month. 

The summer maximum of IASI NH 3 observed in Mari- 
copa has been attributed to the strong dependency of am- 
monia volatilization on surface temperature (Clarisse et al., 
2010b). The similarity of NH 3 and dust cycle suggests also a 
dependency on wind speed on NH 3 volatilization. 

In Africa, there is a distinct peak of IASI NH 3 in May 
which was also observed in surface concentration of NH 3 in 
Banizoumbou by Galy-Lacaux et al. (2001). They related the 
high concentration of ammonia in Banizoumbou to a local 
density gradient of animals grazing in pasture areas, but they 
did not discuss the origin of this pronounced peak in May. 
Except for this peak, the seasonal cycle of NH 3 in Africa is 
quite similar to the one of dust. The highest concentration of 
NH 3 is observed in India. The high level of NH 3 concentra- 
tion in India compared to the rest of the world has already 
been shown by Dentener and Crutzen (1994) and later con- 
firmed by Clarisse et al. (2009). The similarity of dust and 
NH 3 seasonal cycles at all these different sites, except for 
Karachi, suggests a common origin, such that NH 3 may be 
used as a tracer for anthropogenic dust. Next, we would like 
to detect potential differences in optical properties between 
natural and anthropogenic dust, in particular concerning the 
mass extinction efficiency. 


6 Mass extinction efficiency of dust 
6.1 Dependency on NH 3 

Theoretical calculation using Mie theory indicates that the 
mass extinction efficiency € of dust decreases when mixed 
with sulfate or nitrate (Bauer et al., 2007). As we have seen 
that NH 3 is a good tracer for anthropogenic dust, the influ- 
ence of anthropogenic dust on e should be correlated with 
NH 3 concentration. 

To examine the existence of such effect on optical 
properties retrieved from sunphotometers, we extract all 
AERONET inversion Level 2 daily data from April 2009 to 
March 2010. As we are only interested in dusty cases, we re- 
strict the data for days when a is lower than 0.5. To calculate 
the mass extinction efficiency we multiply the volume con- 
centration retrieved from the inversion of AERONET data 
(Dubovik and King, 2000) by the mass density of uncoated 
dry dust (2600 kg m -3 ) and divide by the extinction mea- 
sured directly by the sunphotometer. We have indicated in 
Sect. 3 a possible 10 % error in the density by considering un- 
coated dry dust. We selected the extinction at 670 nm because 
there are more instruments measuring at that wavelength than 
at 550 nm. The fine and coarse mode mass extinction effi- 
ciencies are calculated separately, using the corresponding 
retrieved volume concentration. These values are then clus- 
tered as a function of co-located (same day and within 25 km) 
IASI NH 3 concentration. 

Figure 5 shows the daily values of mass extinction effi- 
ciency 6 of fine (FMEE) and coarse (CMEE) modes at all 
AERONET sites with Angstrom exponent less than 0.5, as a 
function of NH 3 column burden and season. The mean values 
of FMEE and CMEE for different ranges of NH 3 concentra- 
tions are also provided in the Figure. We see a clear decrease 
of CMEE as NH 3 increases, with mean values at 0.62 m 2 g -1 
for NH 3 less than 2mgm -2 and decreasing to 0.48 m 2 g -1 
for NH 3 greater than 20mgm -2 . On the other hand, there 
is a factor 3 variability of FMEE but apparently unrelated to 
NH 3 concentration. Considering NH 3 as a tracer of anthro- 
pogenic dust, these results suggest that the mass extinction 
efficiency of anthropogenic dust has a lower coarse mode 
mass extinction efficiency than natural dust, but has similar 
fine mode mass extinction as natural dust. So, the value of 
€ = 0.6m 2 g _1 derived from Eq. (1) corresponds quite well 
to the mean value observed by AERONET sunphotometer in 
the case of natural dust. 

Figure 5 shows that the variability of FMEE at Zinder Air- 
port in summer is as large as a factor 3, still with the same 
range of NH 3 burden. Similarly, FMEE at Jaipur in spring 
varies by more than a factor 2 while NH 3 is ranging between 
2 and 4mgm -2 . There is clearly another factor controlling 
FMEE variability. 

Sullivan et al. (2007) have shown that ammonium sul- 
fate will preferentially accumulate on fine mode dust be- 
cause H 2 SO 4 reactions with dust depends on aerosol surface 
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Fig. 5. Daily mass extinction efficiency (m 2 g _1 ) for coarse 
(CMEE) and fine (FMEE) modes at all AERONET stations from 
April 2009 to March 2010 when the Angstrom exponent is less 
than 0.5, and as a function of IASI NH3 column burden for 8 
ranges of values (colors), with mean values of FMEE and CMEE 
given in units of m 2 g ~ 1 provided in upper right. The correspond- 
ing season is indicated by star (December- January-February), trian- 
gle (March- April-May), square (June-July-August), and diamond 
(September-October-November) symbols. The Zinder Airport and 
Jaipur AERONET sites are indicated by the letter Z and J, respec- 
tively. 


Mass extinction efficiency and IASI NH 3 
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20 < NH3 < 40 mg m 2 
CMEE=0.48 FMEE=5.66 
10 < NH3 < 20 mg rrr 2 
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8 < NH3 < 10 mg m' 2 
CMEE=0.54 FMEE=6.04 
6 < NH3 < 8 mg rrr 2 
CMEE=0.52 FMEE=6.37 
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1 < NH3 < 2 mg rrr 2 
CMEE=0.62 FMEE=5.35 
0 < NH3 < 1 mg m 2 
CMEE=0.61 FMEE=6.71 



area to which fine mode contributes more than coarse mode. 
They further explained that these reactions will displace any 
preexisting nitrate back to the gas phase and will subse- 
quently accumulate on the coarse mode. Ammonium salts 
were found to be internally mixed with dust, and are strongly 
correlated with the sum of acids in dust. Thus they demon- 
strated that ammonium predominantly accumulates in fine 
mode dust, and subsequently on coarse mode with increasing 
ammonium. As ammonium salts are hygroscopic particles, if 
they are mixed with dust then FMEE should be sensitive to 
relative humidity, which may explain the large variation of 
FMEE at one location and during the same season. 

6.2 Dependency on RH 

Hygroscopic growth of ammonium sulfate can increase e by 
an order of magnitude from its deliquescence point to high 
(> 90 %) relative humidity (Tang, 1996). We should also ex- 
pect some growth of fine dust internally mixed with sulfate, 
although not as pronounced as pure sulfate because dust itself 
is hydrophobic. 
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Fig. 6. Daily mass extinction efficiency (m 2 g _1 ) for coarse 
(CMEE) and fine (FMEE) modes at all AERONET stations from 
April 2009 to March 2010 when the Angstrom exponent is less 
than 0.5, and as a function of relative humidity at 925 mb (upper 
panel) and AERONET imaginary part of refractive index at 670 nm 
(lower panel). The range of RH and k values and the corresponding 
mean values of FMEE and CMEE, in units of m 2 g -1 are provided 
on the upper right of each panels. The Zinder Airport and Jaipur 
AERONET sites are indicated by the letter Z and J, respectively. 


Opposite to natural dust sources, located in arid or semi- 
arid regions, most anthropogenic sources are influenced by 
moist air from monsoon flow (Ginoux et al., 2012), which 
may generate hygroscopic growth of aerosols. But if these 
sources are instead influenced by dry winds or dust is in ele- 
vated layers, no growth should be expected. 

To verify that fine mode dust is mixed with NH3 and 
grows hygroscopically, we extract the daily boundary layer 
(925 hPa) relative humidity (RH) from NCEP/NCAR re- 
analysis (Kalnay et al., 1996) at the AERONET sites for the 
corresponding days and locations. Figure 6 (upper panel) is 
similar to Fig. 5 but each data is color coded as a function 
of RH. The Figure shows no apparent dependency of mean 
FMEE on RH although RH varies below and above the deli- 
quescence point of (NH4)2S04 (Haywood and Ramawsamy, 
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1998). At Zinder airport, the relative humidity in summer 
does not vary considerably ranging from 85 and 90 %. Sim- 
ilarly, at the Jaipur site in spring, RH is below the deliques- 
cence point as RH values varies between 75 and 80 %. These 
results indicate a weak dependency of FMEE on hygroscopic 
growth in these measurements, but it does not mean that there 
is none. 

Several reasons may hinder mesurable changes of FMEE 
by RH. These include the dependency of the deliquescence 
of ammonium salts on hysteresis (Martin et al., 2000, 2008), 
and the presence of other compounds, such as organics 
(Wang et al., 2008), or that dust and NH 3 are in different 
vertical layers. We have also to consider that monsoon con- 
ditions are often accompanied by clouds, which are discarded 
from AERONET Level 2 data by the cloud screening algo- 
rithm. In addition, fine mode mass extinction efficiency is 
mostly sensitive to change in scattering and absorption prop- 
erties rather than size. The opposite is true for the coarse 
mode mass extinction efficiency (Seinfeld and Pandis, 1998). 
So, it may be more effective to study variability of FMEE 
based on absorption property. 

6.3 Dependency on k 

Mass extinction efficiency of aerosol depends on the refrac- 
tive index, whose imaginary part is related to its absorption 
property. Figure 6 (lower panel) shows the dependency of 
FMEE and CMEE on the imaginary part of the refractive in- 
dex (k) at 670 nm. As absorption increases FMEE decreases 
from about 12 to 5m 2 g _1 . The decrease of FMEE with in- 
creasing k appears clearly at the two AERONET sites (Zinder 
airport and Jaipur). 

Analysis of dust samples collected at different locations 
in North Africa indicates k values at 655 nm ranging from 
0.0016/ to 0.0065/ (Wagner et al., 2012), which is cover- 
ing most of our results. The size resolved composition of 
the soil samples by Wagner et al. (2012) indicates a much 
larger abundance of sulfate and other non-crustal elements 
in the fine mode, with sulfate volume abundance varying 
from 2 to 6 %. The authors related the variability of k to 
mineralogical composition. In fact, mineralogy determines 
not only the amount of iron, and thus absorption by dust, 
but also the chemical processing by ammonium salts. Sul- 
livan et al. (2007) have shown that sulfate is more strongly 
associated with aluminosilicate-rich dust particles while ni- 
trate and chloride are more associated with calcite-rich dust. 
Iron is essentially associated with aluminosilicates. This dou- 
ble influence of mineralogy may explain the variability of 
FMEE. The strong dependency of k to source location has 
been shown for Saharan dust using backtrajectory analysis 
(McConnell et al., 2010). Here, we show that such varia- 
tion is also linked to the anthropogenic origin of dust due to 
its mixing with ammonium salts or other acidic components 
from organic soils used for agriculture. 


We can conclude from Figs. 5 and 6 , that over source ar- 
eas, natural dust mass extinction efficiency is 0 . 6 m 2 g _1 for 
coarse mode, and decreases to 0.4-0.5m 2 g _1 for anthro- 
pogenic dust, mostly due to mixing with ammonium salts or 
organics from fertilizer. On the other hand, for fine mode dust 
mineralogy is a major factor controlling the mass extinction 
efliciency. 


7 Conclusions 

The global distributions of dust and NH 3 column burden have 
been compared to determine: ( 1 ) if they are mixed over dust 
sources, ( 2 ) the global and regional fractions of dust mixed 
with NH 3 , and (3) the effect of such mixing on dust optical 
properties. 

To achieve these objectives, dust optical depth (DOD) de- 
rived from M-DB2 aerosol products has been converted into 
column burden by dividing DOD by a theoretical mass ex- 
tinction efficiency. The distribution of M-DB2 dust column 
burden is then compared to NH 3 retrieved from the IASI 
satellite data. 

Using two different land-use datasets we found substantial 
similarity between the annual mean distribution of dust and 
NH 3 burdens, particularly over croplands. The mass of dust 
mixed with NH 3 represents 7 to 8 % of the global mass esti- 
mated to be 4.1 Tg over dust sources, depending on the land- 
use dataset. Over croplands mixed dust represents 22 % glob- 
ally, with only 1 % difference between land-use datasets. On 
a seasonal basis, the maximum mixing is observed in spring 
with almost 30% of dust over cropland mixed with NH 3 . 
Among the different continents, the Indian subcontinent has 
the largest fraction of dust mixed with NH 3 (26 %). Our re- 
sults indicate that NH 3 is a good tracer for anthropogenic 
dust. 

Comparison of dust column burden retrieved from M-DB2 
and sunphotometers at 10 AERONET sites, spread over dif- 
ferent continents, indicates good agreement in Africa and 
Asia but large overestimation of M-DB2 values over Arizona 
and to a lesser extend in Israel. 

To study the difference between anthropogenic and natural 
dust on optical properties, we derived the mass extinction ef- 
ficiency for the coarse and fine mode at all AERONET sites. 
Only dusty days with Angstrom exponent less than 0.5 were 
considered. By collocating IASI NH 3 column burden with 
AERONET data, we found that € is 0.6 m 2 g -1 for the coarse 
mode, and it decreases to a minimum of 0.4-0. 5 m 2 g -1 when 
mixed with concentration of NH 3 in excess of 20mgm -2 . 
This represents a 20 % decrease with a ten-fold increase of 
NH 3 concentration. The mass extinction efficiency for fine 
mode is ten times higher than for coarse mode of unmixed 
dust, with no clear dependency on NH 3 concentration or rel- 
ative humidity. Instead, we found a sensitivity to the imag- 
inary part of the refractive index that we attributed to min- 
eralogy and mixing with other aerosols. In view of previous 
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analysis of in-situ data, we related this decrease of coarse 
mode € as a result of the accumulation of purely scattering 
ammonium, while for the fine mode the amount of absorp- 
tion is determined by the mineralogical composition. Miner- 
alogical composition determines not only the amount of ab- 
sorptive iron but also the accumulation of ammonium salts. 

The relatively weak change of mass extinction efficiency 
is in agreement with a previous model study (Bauer et al., 
2007). However, such study considered dust aging during 
transport whereas we have considered here dust co-located 
mostly over croplands with IASI NH 3 hot spots. The exis- 
tence of mixed dust with NH 3 over large areas implies that a 
significant amount of dust is already mixed with ammonium 
salt before its long range transport and cloud processing. This 
in turn will affect dust lifetime, and its interactions with ra- 
diation, cloud properties, and ocean biogeochemistry. 
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